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Introduction 

Fractional processes are one of the main tools for modeling the phenomena 
of long-range dependence in natural sciences, financial mathematics, telecom- 
munication networks etc. Due to the role played by Gaussian distribution, 
the most popular and the most intensively investigated fractional process is 
the fractional Brownian motion , a centered Gaussian process with the 
covariance function E ^B^ (t)B^ (s) ] = ^(t'^^ + s"^^ ~ 1^ ~ "^l^^)- The param- 
eter H G (0, 1) is called the Hurst parameter and measures the smoothness 
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of trajectories of the process (it is approximately the Holder exponent of the 
process) and the "depth of memory" of the process (for H > 1/2 the process 
exhibits the property of long-range dependence). 

From the point of view of possible applications, there are two main draw- 
backs of fractional Brownian motion. The first one comes from the Gaussian 
distribution, which has extremely light tails, though many data coming from 
applications are heavy-tailed. The second one is the homogeneity of incre- 
ments that does not allow to model processes having different regularity and 
different time dependence properties at different time instances. A related 
problem is a self-similarity property, which briefly means that the proper- 
ties of the process are the same under each scale. However, the absence of 
such property is apparent in many cases and mostly evident in stock price 
processes: long-term data is much smoother than wild intraday quotes. 

The light tails problem is worked around usually by considering fractional 
stable processes. In contrast to the Gaussian case, where the covariance struc- 
ture determines the whole distribution of a process, so there is essentially one 
fractional process, in stable case there are many of them: linear fractional sta- 
ble process, harmonizable fractional stable process, Liouville stable process 
etc. (See book [15|] for an extensive review of different fractional processes.) 

In turn, the homogeneity problem is solved by considering mult if r actional 
processes. Recently, several multifractional extensions of fractional Brown- 
ian motion were defined, based on different representations of the fractional 
Brownian motion: moving average (linear) multifractional Brownian motion 
[l3|], Volterra multifractional Brownian motion 14 1, harmonizable multifrac- 
tional Brownian motion pj]. 

In this paper, we consider a process called real harmonizable multifrac- 
tional stable process which has both properties of heavy tails and multifrac- 
tionality, which can be regarded both as a multifractional generalization of 
a harmonizable fractional stable process and as a stable generalization of 
harmonizable multifractional Brownian motion, and can be used to improve 
models involving either kind of processes. 

Our main interest in this paper is in path properties of this process: con- 
tinuity, existence and joint continuity of local times. For fractional harmo- 
nizable stable process continuity was proved in j?! and local times properties 



were considered in [16| . For related results with multifractional harmonizable 
Levy processes we refer to j^. 

The paper is organized as follows. In Section 1 we give necessary pre- 
requisites on stable distributions and local times. Section 2 focuses on path 
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properties of the process considered: almost sure continuity and localizability. 
Section 3 is devoted to existence and properties of local times. 

1. Pre-requisites 

1.1. Stable random variables and processes 

In this paper we focus only on symmetric a-stable (SaS) random variables 
with a e (1,2). We recall that a random variable ^ is called SaS with a scale 
parameter cr"^ if it has a characteristic function 

An important tool to construct stable random variables is independently 
scattered rotationally invariant complex SaS random measure with the Lebesgue 
control measure, which is a complex-valued cr-additive random measure M = 
Ma on M defined by the following properties. 

1. {Rotationally invariant complex SaS) for any Borel set A C M and 
any G R the distribution of e^^M{A) is the same as of M{A), and 
Y{bM{A) is SaS with the scale parameter \{A). 

2. {Independently scattered) for any disjoint Borcl sets Ai, . . . ,An C [0, oo) 
the values M{Ai), . . . , M{An) are independent. 

3. For any Borel set A C M M{-A) = M{A). 

For a function / : R ^ C such that 

f{-x) = J{x) for all X e R 

and 

= / \f{x)rdx < oo 
■Jr 

it is possible to define a stochastic integral 

/ f{x)M{dx), 
Jr 

which appears to be a real SaS random variable with the scale parameter 

II/IIl"(k)- 

In other words, stochastic integral gives an isometry between the space 
of SaS real random variables spanned by the measure M with the norm 

llell^ = c(0 = -logE[e^«] 



1.1 
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and the subspace of L"(R) consisting of functions satisfying fll.ip . i.e. having 
adjoint values at symmetric points. 

We end this subsection with the so-called LePage representation of pro- 
cesses given as transformations of SaS random measure. For details see 



Assume we have a measurable function /: x M — )■ C such that for 
each t > the function f(t,-) satisfies (II. ip and belongs to L"(]R). Define a 
process {X(t),t > 0} by 



The next proposition is a slight modification of [8j , the proof is exactly the 
same as there with a slight adjustment for the property M{—dx) = M{dx) 
in our case, so we skip it. 

Theorem 1.1. Let ip be an arbitrary probability density on M equivalent to 
the Lebesgue measure. Also let [Tk^k > 1}, {^k,k > 1}, {gk,k > 1} be three 
independent sets of random variables, such that 

• {rk,k > 1} is a sequence of arrivals of Poisson process with unit in- 



• {iki k>\} is a sequence of independent random variables with density 



• {dki k>l} are independent rotationally invariant complex Gaussian 
with E [ |Re ] = 1. 

Then the process t > 0} defined by (\1.2\\ has the same finite- dimensional 

distributions as the process 



AH- 





tensity; 



oo 





where Ca = (/^ 
each t. 



oo 



X ° sin X dx) 



l/a 



and this series converges almost surely for 
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1.2. Local times 

Let X = {X(t),t > 0) be a real- valued separable random process with 
Borel sample functions. For any Borel set B C the occupation measure 
of X on 5 is defined by 

Hb{A) = A({s e B, X{s) e A}) for all Borel sets A in R, 

where A is the Lebesgue measure on R+. If /xb is absolutely continuous with 
respect to the Lebesgue measure on R, we say that X has a local time on 
B and define its local time, L{B, •), to be the Radon-Nikodym derivative of 
fiB- We write L(t,x) instead of L{[0,t],x) and interpret it as the time spent 
by X in x during the time period [0,t]. 

By standard monotone class arguments we deduce that the local times 
have a measurable modification that satisfies the following occupation density 
formula : for any Borel set B C R+ and any measurable function / : R — >■ R 



f{X{t))dt^ / f{x)L{B,x)dx. 



By applying this formula to f{x) = e*""^, and writing L{B,u) for the 
Fourier transform of L{B,x), we get L{B,u) = J^e'^'^^'Us and by the 
Fourier inversion formula L{B,x) = ^/ J^e^'^^-^^^^'^^dsdu, if this integral 
exists. 

As a consequence, the following expressions for the moments of local time 
hold : for any x,y & R, t, h & R+ and m ^ 2 



E[{L{t + h,x)-L{t,x))'"] 



-i— [ I exp(-ia;f;M,)E 



exp < i UjX{i 



JJ^ d,U j dsj^ 
3=1 i=i 



and for every even m ^ 2 

E [ {Lit + h,y)- Lit, y) - Lit + h,x)+ Lit, x))™ ] 



(27r) 



[t,t+h] 



m I mm 



exp<{i^WjX(s^- 



i=i i=i 



Suitable upper bounds for these moments imply the existence of a (jointly) 
continuous version of local time and, as a consequence, a certain degree 
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of irregularity of the sample paths of the process itself. In order to prove 
the joint continuity of the local time of Gaussian processes, S.M. Berman 
(see e.g. jsf) has introduced the notion of local nondeterminism (LND). 



This notion has been extended to stable processes by J. P. Nolan 12|, where 
the equivalent notion of locally approximately independent increments was 
introduced. These notions will be recalled in section 3. Since then the local 
time of stable processes has been studied by several authors; we refer to 



the recent survey by Y. Xiao [l7j for more recent results. The local time of 



multifractional moving average stable processes has been studied in [sj. 

2. Definition and pathwise properties of real harmonizable stable 
process 

Let M be an independently scattered rotationally invariant complex SaS 
measure on M defined in Subsection II. 1[ 

Throughout the paper we will denote by C any constant, which does not 
depend on any variables, unless otherwise is stated. Of course, C may change 
from line to line. 

Recall that a real harmonizable fractional stable process with Hurst pa- 
rameter H is defined as 



Z'^it) = [ ^^^^Midx). (2.1) 

JR \x\ 



A multifractional generalization of this definition consists, naturally, in let- 
ting the Hurst parameter depend on t. 

Definition 2.1. A real harmonizable multifractional stable process (rhmsp) 
with Hurst function H(t) and a stability parameter a is defined as 

m = I -^^^M{dx). (2.2) 

Clearly, X{t) = Z^^'\t). We assume that < H = inU H{t) < H{t) < 
snpt H{t) = H <1. 

2.1. Norm estimates for the increments 
Lemma 2.2. For all Hi,H2 E {H,H) it holds 

ll^^^i(t) _ z^^(t)|| <c\H^-H2\, te [o,r]. 
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Proof. Write 

^C\Hi-H2r([ \x\''^^-^^-^\\og\x\\°'dx+ [ \x\-^-°'^\\og\x\\''dx) 

\J\x\<l J\x\>l / 

<C\H,-H2\'', 

whence we have the assertion. □ 

We assume that H is Holder continuous with order greater than H, i.e., 
there exists ^ > H s.t. for all s > 

\H{t)-H{s)\ < C\t-s\'^. 

Lemma 2.3. There exist positive constants Ci,C2 > such that for any 
He [H,H] one has 

Ci\t - < ||Z^(i) - Z''{s)\\^ < C2\t - s\" 

locally uniformly in s,t. 

Proof. Write 

\\Z^{t) - Z"{s)\\l = I |e**^ - e^^^P \x\-^-''"dx <C [{lA\t-s\ \x\y\x\-^-'^^dx 

V J\x\<l/\t-s\ J\x\>l/\t-s\ J 

<C{\t- sl"-""^^-"^ + \t- sr'' )=C\t- si""" . 

To prove the lower bound, observe that there exist positive constants Ci,C2 
such that |e*^ — 1| > Ci for \y\ < C2 and write 



\z^it)-z"is)\\:> [ 

J \x\ 



'\x\<C2/\t-s\ 

>c\t-sr I ixi^'^^-^^-'dx ^c\t- . 

|a;|<C2/|t-s| 

□ 
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Lemmata 12.21 and 12.31 imply the following 



Proposition 2.4. There exist 6,Ci,C2 > s.t. for rhmsp X given by fl2.2l) 
and \t — s\ < 6 it holds 



Ci\t 



H{t,s) 



r'''^ <\\X{t)-X{s)\l<C2\t-s 
where H(t, s) = min^t s] H{u), H(t, s) = max[f H{u 
Proof. Let H{t, s) = H{t), H{t, s) = H{t). 

\\m'X{s)\L 



t,se[0,T] (2.3) 



< 



ZHis), 



ZH(t,s), 



Z^^'''\t) - Z^^'''\s) 



< \H{t) - H{i) \ + \H{s) - H{i)\ + C\t - < C\t - sp + C\t 

Since H(t, s) < 7, we get the upper bound. 
The lower one is proved similarly: 



C\t 



\H{t,s) 



< 



< \\X{t)-X{s) 



a 



ZH{t,s), 



< \\X{t) - X{s)\l + \H{t) - Hii)\ + \H{s) - H{i)\ 
<\\Xit)-Xis)\l + C\t~s\\ 



□ 



2.2. Holder continuity of rhmsp 

In this subsection we prove a Holder continuity of rhmsp. Our argument 
is a slight modification of the one found in 0] for harmonizable fractional 
stable motion. A similar argument was also used in [^[| to prove a Holder 
regularity of operator scaling stable random fields. 

Theorem 2.5. The rhmsp X has a version, which is almost surely Holder 
continuous of any order n < H and moreover almost surely satisfies 



sup \X{t)-X{s)\=o{6"\\og6 

t,se[o,T] 

\t-s\<S 



l/a+l/2+e^ 



5^0+, 



for all T,e>0. 
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Proof. Let T > be fixed and throughout this proof t, s G [0, T]. 

We use the LePage representation fll.3p . To simphfy the notation we 
write this representation for the process X itself rather than for its version: 

X{t) = C„Re5^r-%(a)-^/"/(t,a)^/c, 
fe>i 

where f{t,x) = (e**^ — 1) \x\^^^"'^^^^\ (p{x) = Kr^\x\~^ |log |a;| 1""'"^'', ?7 > is 
arbitrary but fixed, K^^ is a normahzing constant. 

Conditioning on F and C,, X has the Gaussian distribution, so 

E [(X(t) - X{s)f I r,e] = C^5^r,^/V(efc)-^/° - < Ca{u), 

k>l 

where 

aH = J]r-^/>(6)-2/" sup \fit,Q-fis,^,)f. 

k>l \t-s\<u 

Write 

sup \f{t,^k) - f{s,^k)\ 

\t~s\<u 

< sup |(e^*^ - e^'^l Ixp^/""^^*^ + sup |e^"^ - l| jxp^/" jx^^^*) - 

\t—s\<u \t~s\<u 

< C{u \x\ A 1) \x\~^^" {\x\~^ V \x\-^) 
+C(|x| Al)|,Tr^/"(|xr^V|a;r^)|log|x|| sup \H{t) - H{s)\ 

\t—s\<u 

< C Ixp^/" (Ixp^ V |xr^)((M |x| A 1) + A 1) |log \x\\u"'). 

Keeping this estimate in mind, take now the expectation [a{u) ] with re- 
spect to the variables ^ only: 

E^[a{u)]<CS{T){h + l2), 
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where 

Ji= / \x\-^/''{\x\-^'' y\xr^''){u\x\Mf^^-^/''{x)dx 
' x-\x-^^ V x-^^){ux A If Ilogxl^^+^^^'Z"-^) dx 

Jo 

/•oo 

= 2u^^ / V x-^^){x A 1)2 |logx|(^+'')('/"-^)+' rfx < Cu^\ 

Jo 

S{T) = ^r;'/" < oo a.a. T, 

fc>i 

where the last is true owing to the fact that Tj/j — )• 1, j — )• oo, almost surely 
by the strong law of large numbers, and 2/ a > 1. Therefore 

[a{u) ] < C{ry^ |logn| 

almost surely. 

Define b{u) = u^^ |logMp^+'')/°. We have 



n>l 



6(2-") 



<C{T)Y,n-' 



5 



n>l 



SO for almost all ^, F we have a(2 "')/6(2 '^) — )■ 0, tt, — )■ cxd. It is easy to see 
that b{2t) < Cb(t), and a{u) is increasing, so from the last convergence we 
get a{u)/b{u) — )■ 0, m — )■ 0+, or a{u) = o^^riu^^ jlog-ul^''^'''^'*''"). So we have 

E[{X{t)-X{s)f I r,e] = o^,r{u'^ llognl^^^^")/"), n ^ + . 
Now recall once more that X is Gaussian given ^ and F, so by Lemma 1 of 



sup |A'(«)-X(s)|=o„(i''|log«|"''+''''+"'), «^0+, 

\t~s\<5 

whence we get the statement of the theorem. □ 
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2. 3. Localizahility of rhmsp 

We start this section by giving Falconer's notion of localizahility. 

Definition 2.6. Process X is called H-localizahle at a point t with the local 
version Y if 

|^(X(t + 5M)-X(t)),M>o| ^{r(u),?z>0}, (5^0 + . (2.4) 

(Here stands for the convergence of finite- dimensional distributions.) 

It is called strongly if-localisable at a point t if in (12. 4p the convergence 
is in the sense of the distribution on the path space. 

Some authors use the term local asymptotic self- similarity for localiz- 
ahility, which reflects the fact that the local version Y is an if-self-similar 
process. 

Theorem 2.7. The rhmsp X is localizable at any point t with local version 
being real harmonizable fractional stable process with Hurst parameter H{t). 

Proof. Define 

Yf{u) = ^iXit + 6u)-Xit)). 

We will assume throughout that 6 < 1. 

For ui, . . . ,Un > 0, Ai, . . . , A„ G M denote = t + 6uk and write 



-logE 



exp <^ zV" AfcF/(nfc) 



Now estimate the integrand multiplied by \x\: 



Hit) 



k=l 



k=l 



Hit) 



-Hisk) _ ur-f^W 



-aHit) \^iskX _ ^itx I 



'jb\ 

fe=i 

< C E lA^r (|log \x\r kl""' - H{tX (1 A + |a;r"^(*) js^ - tj" (1 A 

fe=i 

n 

<CY, |A,r {\xr" V A [ - tr |log |a;|r + - C] 

fe=i 

n 



k=l 



< \xr^'-"^ i|,i<i + i^r"^ i|,i>i) (1 + |iog 



x\ 



We remark that the constants here depend only on a, t,ui, . . . ,Un and 
Ai, . . . , A„. 
Now 



1 1 



gaHit) 



X\ 



n 
k=l 



itx 



\X 



<C{\ 



X 



ail-H)-l 



I|a:|<l + \x\ 



Hisk) 



-l-aH 



\X 



Hit) 



I|^|>l)(l + |log 



x\ 
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which is integrable over R. Hence by the dominated convergence theorem 



hm log E 





i 1 


/ -. — r hm < 





cxp <{ i ^ XkY^\uk) 

k=l 



k=l \-^\ 



itx 



X 



itx 



\X 



H{t+Suk) 



\X 



r 1 _ _ rpiukV _ „-ity/S 



dx 



d{5x) 



k=l 
1 



EAfc hm 



fe=i 



\y 



H{t) 



dy 



dy, 



where 
Rs — — 

\y 

Estimate 



1 _^H{t+5uk)-H{t) _ 1 



H{t+6uk) 



\H{t) 



S ^H{t+5uk)-H{t) 



- 1 



log ( n 



\y\ 



< \H{t + 5uk) - H{t)\(\\og5\ + \\og\y\ 
<CSmiogS\ + \log\y\\) ^0, 6^0 + . 



Thus Rs^0,5^ 0+. 
FinaUy, 



hm log E 

5-^0+ 



exp <i^XkY^\uk) 



k=l 
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-J 


R \y\ 



fe=i 



ginfe?/ _ 



dy, 
(2.5) 



which is exactly the logarithm of the characteristic function of . . . ,Z^^*\un), 

as required. □ 

Remark 2.8. By using the same kind of argument as the one used in the 
proof of continuity, it is possible to prove tightness of laws of processes on 
the space of continuous paths and whence derive a strong localizability. 
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3. Local times for rhmsp 

3.1. Properties of the local time 

We start this section by showing the existence and square integrability of 
a local time. 

Proposition 3.1. The rhmsp X has a square integrable local time L{t,x). 

Proof. According to [5], it is enough to check the following "condition (H)": 

s\ < p 



there exists p > and H G (0, 1) and G L^(]R) such that for all \t - 
|E [exp{iX{X{t)-X{s))}] \ <^(A|t-s|^). 
But, in view of (12.31) . for t and s close enough 

|E [exp {^X{X{t) - X{s)) } ] \ = exp {- \Xr \\X{t) - X{s) fj 



(3.1) 



<exp<^-C|A|°||t 



<exp<^-C|A|"||t 



□ 



whence we have (13. ip with = exp {— C|2;|"}, H = H. 
In order to prove further properties, we need 

Definition 3.2 (jlij). A stable random process X is ||-||^ locally non- deterministic 
(LND) on T if 

(LI) ||X(t)||^ > for all t e T; 

(L2) ||X(t) - X{s)\\^ > for all sufficiently close distinct s,t G T; 
(L3) for any n > 1 there exists C„ > s.t. for any ti < t2 < ■ ■ ■ < tn E T 
sufficiently close together one has 

||X(t„) - span {X(ti), . . . , X(t„„i)} t > C„ ||X(t„) - X(t„_i) ^ • 

(3.2) 

In 12| it is shown that the local non-determinism property is equivalent 
to the property of ||-||^ locally approximately independent increments, which 
consists of properties (LI), (L2) above and 

(L3a) for any n > 1 there exists C„ s.t. for any ti < t2 < ■ ■ 
sufficiently close together and any oi, . . . , a„ G M one has 

n-1 

aiX(ti) + ^ afe (X(4+i) - X{tk)) 



< t„ e T 



fc=i 



n-1 



(3.3) 



> a i|aiX(ti)iL + E \Mx{t,+i) - x(tfc)) 11^ 



k=l 
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Theorem 3.3. For any e > the rhmsp X is LND on [e,T]. 

Proof. The main difficuhy is to prove property (L3) of LND, as property 
(LI) is obvious and property (L2) follows from (12.31) . 
We proceed in two steps. 

Step 1. We prove LND for a modification of rhmsp X defined by 



Y{t) 



1 -e~**^)(-ia;) 



-H{t)~l/a 



M{dx), 



(3.4) 



where 



\X\ 



-K giTrK signa;/2 



The Fourier transform of the function fY{t,x) 
(w.r.t. the second variable) on L"(]R) is 

1 



{t - x) 



H{t)-1/I3 



;i - e-^*^)(-2x)-^w-i/" 
-xrr ' I , (3.5) 



T{H{t) + 1/a) 

where P = a/{a — 1) is the exponent adjoint to a, see Lemma [A. 2[ 

In order to check property (L3) for Y, we have to find a good lower bound 

to 



n-l 



Y{t^)-J2^kY{tk) 



k=l 



n-l 



fritn, ■) — ^ Ukfyitk, ■] 



k=l 



It is fortunately given by the Hausdorff- Young inequality: 



n-l 



k=l 



> C 



fritn, ■) — ^'Wfc/y(tfc, 
fvitn 



> C 



n-l 



fritn, ■) — y^^Ukfritk, 



k=l 



c 



V{H{tn) + 1/a) 



{tn - X)'3^(*")-Mx 



> C(t„ - t„_i)^(*") > C ||X(4) - X(t„_i) 



for tn and tn-i close enough. (We have used the fact that fyitk.x) vanishes 
on for < n in the middle, and inequality (12.31) in the last step.) 

But it is straightforward to check (see a much stronger statement below in 
the Step 2) that 

\\X{tn) - X{tn-l)\\^ > \\Y{tn) - Y{tn-l)L " C\H{tn) " H{tn-l)\ 

> ||F(t„) - r(tn-i)IL - c\tn - tn-iT > c \\Y (tr,) - r(t„-i)IL , 
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which gives the desired LND property. 

Step 2. Here we show how the property of locally asymptotically inde- 
pendent increments for Y implies that for X. Denote fx{t,x) = (e**^ — 
1) |a;r^(*^"^/", fi{t,x) = -frit, -x) = /^(t, a;)e-^'^(^W+i/")«ig"^=^/2 and write 
for < < ^2 < • • • < *n < ^ and ai, . . . , an G M (we put to = for the sake 
of simplicity) 



n— 1 



k=l 



n-1 



^ak{fx{tk+i, •) - fx{tk, •)) 



k=0 



n-1 



-iTTiH{ti)+l/a)si: 



fe=0 



> 



rt— 1 



^ak{fi{tk+i, •) - fi{tk, •)) 



fe=0 



n-1 



k=0 



where 



_ ^-in sign x/ (2a) ^^-i-KH(ti) sign x/2 _ ^-i-w H (t) sign x / j ^-j 



Estimate 



^ ||(g-^.f^(t.)sign./2 _ ^-i.Hit,)signx/2^^j^^^^^^^ .) _ ^^(^^^ . 

< C\H{tk+l) - H{tk)\ \\fx{tk+l,-)\\L-'{R) 

+C \H{tk+i) - H{ti)\ \\fx{tk+i, ■) - fx{tk 



= o{\\X{tk+i)-X{tk)\\J, ^oo. 
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Further, 



n-1 



flitk,-)) 



k=0 



n-1 



^ak{fY{tk 



fvitk,-)) 



n-1 



5^afc(r(tfc+i)-y(tfc)) 



A:=0 



k=0 

n-1 



>CY,hk{Y{h+,)-Y{t,))\l 



k=0 



where the last inequahty is true thanks to LND property of Y. Similarly to 
the first estimate of Step 2, 



- Yit,)\l = ||/y(tfc+i, •) - fvih, OIL^d,) = ||/i(tfc+i, •) 

> ||X(4+i) - Xitk)\l - II A(4+i, ■) - A(4, OIL.pR) 
> ||X(4+i)-X(4)L-o(||X(4+i)-X(4)||J, 

so finally 

n-1 

^afc(X(4+i) -X(4)) 



/i(4,-)IL 



k=l 



n-1 



> 



J2 Wk\ ( ||X(tfc+0 - X(t,)L - o{ ||X(4+i) - X(4) t ) 

n-1 

>cY,M\\{x{tk+i)-xit,))\l 



k=l 



for |t„ — ti| small enough. 



□ 



Remark 3.4. It is possible to make the presented proof shorter by skipping 
several lines at the end of Step 1: in fact the interim lower estimate by 
||X(t„) — X(t„_i)||^ is exactly what is needed in proof, and one does not 
need to go further obtaining the LND for Y. Nevertheless, these lines makes 
the proof more structured, and we think that this intermediate result of LND 
for Y is interesting on its own. 

Thanks to 12|, Theorem 4.1] and estimates for the norms of increments 
of rhmsp X we have whence the following result. 

Theorem 3.5. The local time L{t, x) of the rhmsp X is jointly continuous in 
{t,x) fort > 0, moreover, for any k, < {1/H — 1)/2 it is K-Holder continuous 
in X. 
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Remark 3.6 (Gaussian case). It is not hard to see that all results of this paper 
extend to the Gaussian case, viz a = 2. In fact all the proofs work in the 
Gaussian case as well, except the one of the pathwise continuity, but there one 
has a much shorter and direct proof. Also, in the proof of Theorem 13.31 one 
does not need the Hausdorff- Young inequality, but just Parseval's identity. 

It should be noted that some papers mistakenly claim that in Gaussian 
case, the real harmonizable multifractional Brownian motion is up to a mul- 
tiplicative constant equivalent to moving average (or linear) multifractional 
Brownian motion, as defined e.g. in jsj. The difference between the two defi- 
nitions is shown in jcj to be essential. Nevertheless, this claim is true for the 
process Y defined in the proof of Theorem 13.31 



Appendix A. Fourier transform 

In this appendix we compute the Fourier transform which is used by many 
authors, however, we were not able to find a rigorous derivation. Below we 
define the Fourier transform 

fiu) = [ e^"^/(x)da;, 
Jr 

and use the notation vo. 

Lemma A.l. For h G (1, 2), t > the Fourier transform of 

f^ ^{x) = (1 - e-^*^)(-zx)-^ = (1 - e"**^) e'^^'''^'''''^ 

is 



2tx 



hM = ^Mt 



Proof. First note that = e~'"^v'^~^dv is analytic for Re^; > 0. So 

it follows from the operational calculus that for any a > 

„,h~l 1 na+ico -i poo— ai 

-t- = / e^'^z-^dz = / e-'''n-iy)-^dy, 

T{h) 2m J^_i^ 2tt J _^_a^ 

where we have changed the variable z — )• iy. Plugging v = —u and v = t — u 
to this identity, we get 

-| -I POO— ai 

■{{t - u)'l-' - {-u)1-') = - / e™ni - e-^'y){-zyr'dy, (A.l) 



Tih)"-' V ^+ ^ 27r 
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Now let in this integral y = x — ai, x E W and estimate for a G (0, 1) the 
integrand as 

|e»s/(l _e-%)(_^^)-/^| = ga« 1^ _g-ai-ite| ^y^-h ^ c{u){t\x~ai\ A 1) \yf^ 
= C{u){t \y\-''^^ A lyl-"") < C{u){t \x\-''^^ A Ix^''), 

(A.2) 

which is integrable due to the assumption h G (1,2). So letting a — )■ 0+ in 
f lA.ip yields the desired result by the dominated convergence theorem. □ 

By the Hausdorff- Young inequality (see 0, Theorem 5.7]), for a G [1, 2] 
the Fourier transform from L^(M) fl L°(R) can be extended to a bounded 
linear operator J-q: L"(R) L^(R), where P = a/{a — 1) is the exponent 
adjoint to a. We will call this map a Fourier transform on L°(R), and we 
emphasize once more its boundedness due to the Hausdorff- Young inequality: 

ll-^a/llL'3(R) — II/IIl"(K) • (^-3) 

The following lemma is an L"(R) analogue of Lemma [A.ll 

Lemma A.2. For a G (1,2), h G (1/a, 1 + 1/a) and t > the Fourier 
transform on of 

/^ ,(x) = (1 - e~'''')i-ix)-^ = (1 - e"'*") |a;|-'^e*-'*^ig--/2 

is 

Proof. Repeat the proof of the previous lemma to inequality ( 1A.2P and raise 
it to the power a: 

|e™^(l - e-'*^)(-zy)-'^p < C{uY{t \xf^-^^ A |x|-"'^), 

which is integrable for /i G (1/a, 1 + 1/a). So the integrands in ( lA.ll) converge 
as a — )• 0+ in L"(R) to fh,t{x) by the dominated convergence theorem, hence 
by continuity of J-'q, on L"(M) we get the statement of the lemma. □ 
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